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Resonant annihilation of extremely high-energy cosmic neutrinos on big-bang relic anti-neutrinos 
(and vice versa) into Z-bosons leads to sizable absorption dips in the neutrino flux to be observed 
at Earth. The high-energy edges of these dips are fixed, via the resonance energies, by the neutrino 
masses alone. Their depths are determined by the cosmic neutrino background density, by the 
cosmological parameters determining the expansion rate of the universe, and by the large redshift 
history of the cosmic neutrino sources. We investigate the possibility of determining the existence of 
the cosmic neutrino background within the next decade from a measurement of these absorption dips 
in the neutrino flux. As a by-product, we study the prospects to infer the absolute neutrino mass 
scale. We find that, with the presently planned neutrino detectors (ANITA, Auger, EUSO, OWL, 
RICE, and SalSA) operating in the relevant energy regime above 10 21 eV, relic neutrino absorption 
spectroscopy becomes a realistic possibility. It requires, however, the existence of extremely powerful 
neutrino sources, which should be opaque to nucleons and high-energy photons to evade present 
constraints. Furthermore, the neutrino mass spectrum must be quasi-degenerate to optimize the 
dip, which implies m v > 0.1 eV for the lightest neutrino. With a second generation of neutrino 
detectors, these demanding requirements can be relaxed considerably. 

PACS numbers: 14.60.Pq, 98.70.Sa, 95.85.Ry, 95.35.+d 



I. INTRODUCTION 

Neutrinos are the elementary particles with the weak- 
est known interactions. Correspondingly, they can prop- 
agate to us through the cosmic microwave and neutrino 
background (CMB and CjaB, respectively) without sig- 
nificant energy loss even from cosmological distances. 
A possible exception to this transparency is resonant 
annihilation of extremely high energy cosmic neutrinos 
(EHECi/s) on big-bang relic anti-neutrinos (and vice 
versa) into Z-bosons [1-4]. This occurs near the respec- 
tive resonance energies, 

E™ = = 4.2 x 10 22 eV ( ^1) . (i) 

2m Ui \ m Vi J 

with mz =91.2 GeV denoting the mass of the Z-boson [5] 
and 

(J 1 — I? 2, 3) the non-zero neutrino masses — 
for which there is rather convincing evidence inferred 
from the apparent observation of neutrino oscillations [5] . 
On resonance, the corresponding cross-sections are en- 
hanced by several orders of magnitude. This leads to a 
few percent probability of annihilation within the Hub- 
ble radius of the universe, even if one neglects further 
enhancing effects due to cosmic evolution. Indeed, it 
appears that - apart from the indirect evidence to be 
gained from cosmology, e.g., big-bang nucleosynthesis 
and large-scale structure formation - this annihilation 
mechanism is the unique process [82] having sensitivity 



to the Cid3 [1]. Moreover, observation of the absorp- 
tion dips would present one of the few opportunities to 
determine absolute neutrino masses [9]. However, the 
mechanism requires that there exists a sufficiently large 
EHECi/ flux at the resonant energies of Eq. (1). One of 
the purposes of this work is to quantify how large this 
EHECf flux must be. 

Since the original proposal in 1982 [1], significant ad- 
vances have been made in theoretical and observational 
cosmology, experimental neutrino physics, and EHECf 
physics. Each of these three areas impacts immediately 
on the CiaB measurement. Thus, it is timely to investi- 
gate again the possibility of determining the existence of 
the GvB and to study the prospects for determining the 
neutrino masses via resonant neutrino absorption. 

What are the new findings that affect the EHEC^ ab- 
sorption probability? 

• The original calculation of neutrino absorption was 
done for a matter-dominated flat universe with- 
out a cosmological constant. Recent observations 
of large-scale gravity, deep-field galaxy counts, and 
Type la supernovae favor a universe with energy 
fractions fl\ « 0.7 in the cosmological constant 
and f2jvf ~ 0.3 in (mainly cold and dark) mat- 
ter [10, 11]. The position of the first Doppler 
peak in recent CMB measurements strongly sug- 
gests that the universe is flat, i.e. the fractional 
contribution of curvature energy fi/. is negligibly 



2 



small. These cosmological parameters, together 
with the Hubble parameter H , determine the ex- 
pansion rate of the universe as a function of look- 
back distance. This expansion history, in turn, cru- 
cially impacts the EHECV fluxes at Earth, since 
their sources are almost certainly located at cos- 
mological distances [83]. 

• The oscillation interpretation of the atmospheric 
neutrino data offers a lower limit on the heaviest of 
the three mass-eigenstates of 



in,, 



> 



2 

atm 



> 0.04 eV 



(2) 



at 95% confidence level (CL), from the inferred 
mass splitting Am 2 tm [84]. On the other hand, 
studies of the cosmic evolution of the large-scale 
structure, as observed today, from primordial den- 
sity perturbations, measured in the CMB, yield an 
upper bound on the sum of neutrino masses [10-12] 



<1.2 eV. 



(3) 



Since oscillation studies also reveal that neutrino 
mass-splittings are small compared to the eV 
scale, the cosmic bound per mass-state is conser- 
vatively [85] ~ 0.4 eV - about a factor of five 
better than the laboratory bounds inferred from 
tritium beta decay [14, 15], m U3 < 2.2 eV (95% 
CL) [86], and neutrinoless double beta decay [16], 



< (0.66 4- 2.70) eV [87]. Thus we have 
0.04 eV < m v , < 0.4 eV . 



(4) 



It is remarkable that the neutrino mass, whose 
value was compatible with zero at the time of the 
original proposal for Z-dips [1] , is now known to be 
not only non-zero, but to lie within a one-order of 
magnitude range! Accordingly, the resonant anni- 
hilation energy for the heaviest mass-state is also 
known to within an order of magnitude: 



1 x 10 22 eV < Elf < 1 x 10 



2:S 



eV. 



(5) 



This resonant energy, when divided by the Z-decay 
multiplicity of ~ 40, predicts secondary cosmic ray 
particles with energies spanning a decade or more 
above the Greisen-Zatsepin-Kuzmin (GZK) energy, 
Eqzk = 4 x 10 19 eV. This is the energy beyond 
which the CMB is absorbing to nucleons, due to 
resonant photopion production [19]. The associa- 
tion of Z-bursts with the mysterious cosmic rays 
observed above Eqzk is a controversial possibil- 
ity [20-25]. Intriguingly, the neutrino mass window 
required in such a scenario coincides quantitatively 
with Eq. (4) [24]. 

• Recent proposals for progressively larger EHECf 
detectors, such as the Pierre Auger Observa- 
tory [26], IceCube [27], ANITA [28], EUSO [29], 
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FIG. 1: Current status and next decade prospects for EHEC^ 
physics, expressed in terms of diffuse neutrino fluxes per fla- 
vor, F Va + Fv a , a = e,fi,r; full mixing among the flavors en 
route to Earth [33] is assumed. 

Top: Upper limits from RICE [37], GLUE [38], FORTE [39], 
and Fly's Eye and AGASA [40]. Also shown are projected 
sensitivities of Auger in i/ e , modes and in v T mode (bot- 
tom swath) [41], ANITA [42], EUSO [43], and SalSA [44], 
corresponding to one event per energy decade and indicated 
duration. 

Bottom: Roadmap for improvement in the next decade (2008 
and 2013; adapted from Ref. [32]), corresponding to one event 
per energy decade, as well as the current (2003) observational 
upper bound (solid-shaded) obtained from Fig. 1 (top). For 
the year 2008 (long-dashed), we assume 3 yr of Auger data 
and one 15 d ANITA flight. For 2013 (dashed-dotted), we 
assume 8/3/4 yr Auger/EUSO/SalSA, and 3 ANITA flights. 
The sensitivity will improve if further projects such as Auger 
North and OWL [30] are realized, or if the EUSO flight time 
is extended. Also shown is a wide sample of predictions for 
EHECf fluxes (discussed in § II B). 



OWL [30], and SalSA [31] offer credible hope that 
the collection of an event sample above 10 21 eV 
may be realized within this decade [32] (cf. Fig. 1). 
Another encouraging sign for the progress in ex- 
perimental sensitivity is that existing EHEC^ "ob- 
servatories" , such as RICE [34], GLUE [35], and 
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FORTE [36], have recently put the first sensible 
upper limits on the EHECV flux in the region (5) 
relevant for neutrino absorption [37-39]. We dis- 
play these limits in Fig. 1 (top). 

The organization of this paper is as follows. In § II, we 
determine the EHECV spectra, with their all-important 
absorption features, to be observed at Earth. In par- 
ticular, the depths, widths, and locations of the relic 
neutrino absorption dips are calculated, for various pro- 
posed sources of a diffuse EHECz^ flux. The experimental 
prospects to detect the absorption dips within the next 
decade or beyond are discussed in § III. Finally, in § IV, 
we summarize our results and conclude. 



II. EHECi/ SPECTRAL DIPS 

Given an EHECV source emissivity distribution, one 
can determine the resulting neutrino spectrum to be ob- 
served at Earth, by taking into account resonant annihi- 
lation with the C^B, and energy losses due to redshift. 
We formulate this problem, in § II A, in terms of prop- 
agation functions [45, 46]. We calculate these functions 
by means of the procedure outlined in the original papers 
on relic neutrino absorption spectroscopy [1, 2], updated 
to modern cosmological parameters. The observable neu- 
trino flux arriving at Earth is then obtained by folding 
the propagation function with the EHECz^ source emis- 
sivity distribution. The latter basic input is not known 
in the energy region of interest. Therefore, we intro- 
duce in § II B various parameterizations to model neu- 
trino emission from the most relevant classes of possible 
sources - astrophysical accelerators (bottom-up "Zcva- 
trons" ) or exotic massive particles and topological defects 
(top-down). In § II C, we study the location, depths, and 
widths of relic neutrino absorption dips in the context of 
these model classes of sources, and for various neutrino 
mass scenarios. 



expressed as 



F Va (E) 



d 4 A, y 



dEdAdtdfl 



OO OO 

I /■,„/", -dP Va \ V0 {E;Ei,r) 



4tt 
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dE 



(6) 



C V/3 (r,Ei) , 



where the "propagation distance" r = c t is related to the 
redshift z by 



dz= (l + z)H(z)dr, 
with the evolving Hubble parameter given by 



(7) 



h 2 (z) = Hi [n M (i + zf + n k (i + z f + n A ] . (8) 

Here, H = h 100 km/s/Mpc, with h = (0.71 ± 
0.07) x o!gi [5]> denotes the present value of the Hubble 
parameter. In Eq. (8), 0,m, ^fc, and are the present 
matter, curvature, and vacuum energy densities, respec- 
tively, in terms of the critical density. From the Fried- 
mann equation comes the constraint that fractional en- 
ergies must sum to 100 %, i.e. Qm + ^fe + = 1- As de- 
fault values we choose fi_M = 0.3, Cl k = 0, and SIa = 0.7. 
These values are collectively known as the "concordance" 
model, for they fit a wide assortment of data [10, 11] [89]. 

Justified approximations are that (i) the only type of 
energy loss is due to the redshift, and (ii) the only rele- 
vant interaction is absorption on the relic neutrino back- 
ground, dominated by resonant Z-production [90]. With 
these approximations, the differential propagation func- 
tion is given by (see the Appendix for a thorough deriva- 
tion, properly taking into account neutrino mixing ef- 
fects) 



dP Va \ v AE;Ei,z) 



dE 



(9) 



5 E 



E, 



1 + z 



Y,\U aj \ l P Vj {Euz)\Up j \\ 



A. Neutrino propagation functions 

Let £ Ul3 (r 7 Ei) be the EHECz^ source emissivity dis- 
tribution, i.e. the number of neutrinos vp of flavor fj = 
e,/i, r, per co-moving volume, per unit of time and per 
unit of energy as measured at the source, injected in the 
CVB at a "propagation distance" r = ct from Earth [88] 
with an energy The propagation through the CVB 
can be described [45, 46] by the functions Pb\ a (E; Ei,r), 
which are defined as the expected number of particles 
of type b above the threshold energy E if one particle 
of type a started at a distance r with energy Ei. With 
the help of these propagation functions, the differential 
flux of neutrinos of flavor a (b = v a ) at Earth, i.e. their 
number N Va arriving at Earth with energy E per units 
of energy, area (A), time (t) and solid angle (Q), can be 



where U a j is the leptonic mixing matrix and P V} {Ei, z) is 
the survival probability of a cosmic neutrino i>j injected 
at a redshift z with energy Ei = E (1 + z), 



P Vj (E(l + z),z) = 



(10) 



exp 



dz 



H(z)(l + z) 



with 



s = 2m Vj E(l + z). 



(11) 



We remind the careful reader that z is the emission red- 
shift of the cosmic source, while z is the redshift at the 
time of neutrino annihilation, the latter being integrated 
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from present time back to the emission time. The expo- 
nent in the brackets is (minus) the optical depth (also 
called the "opacity") for a neutrino emitted at redshift 
z. 

The momentum of a massless neutrino today is of or- 
der of the CMB energy, <~ 3 T 7 ~ 0.7 meV. The neutrino 
momentum at an earlier epoch is then ~ 0.7 (1 + z) meV. 
Thus, relic neutrinos will be non-relativistic as long as 
m vj ^ {Pvj ) ~ 0.7 (1 + z) meV holds. In most of the 
region in redshift-space which we consider, the relic neu- 
trinos are non-relativistic. This means that their number 
densities are given by [5] 

n Vj (z)=n Pj (z) = (n,)o (1 + z) 3 = ^(n 7 ) (l + 5) 3 

= 56 cm~ 3 (1 + if , (12) 
with {n v )o and (n 7 )o denoting todays (subscript "0") av- 



P Vj (E(l + z),z) ~ exp 



within the region of support indicated, and identically 
one (no absorption) outside the region of support. Here, 
£™ s is the resonant energy in the rest system of the relic 
neutrinos, given in Eq. (1). 

Numerically, the annihilation probability on the CfB, 
neglecting cosmological expansion, is on the few percent 
level, 

Mo KT) = {0 J1/h) x g o % (16) 

Ho 

Taking cosmological expansion into account, the annihi- 
lation probability is enhanced by the redshift-dependent 
ratio in the exponent appearing in the expression (15) 
for the survival probability. The enhancement is easy 
to understand: In the numerator of the ratio, the factor 
(El es /E) 3 = (1 + z) 3 accounts for the higher target den- 
sity of the CVB in physical volume. The denominator ac- 
counts for the time or path length available per unit red- 
shift for annihilation. What is noteworthy here is that, 
at early times, the Q,m term dominates, and so the neu- 
trino optical depth scales as 1/V^m • Thus, the smaller 
is Om, or equivalently, the larger is the greater is 
the neutrino absorption probability from sources with 
1 + z > (n A /n M ) 1/3 - For example, the ACDM universe, 



erage number density in relic neutrinos and relic photons, 
respectively [91]. 

The annihilation cross-section cr* n ? (s), which is dom- 
inated by the Z-production peak with a rather narrow- 
width, can be approximated by [1] 

H<" n > <5 - 1) , (13) 

where 

/ds 
— = 2ttV2 G f = 40.4 nb (14) 

m z 

is the energy averaged cross-section [92] Gf = 1.2 x 
10~ 5 GcV~ 2 being the Fermi coupling constant [5]. 
Therefore, the survival probability (10), at the injection 
energy Ei = E (1 + z), can be approximated by [1] 



with Q,m — 0.3 and Oa = 0.7, produces absorption 
dips nearly twice as deep as a pure CDM universe with 
Q.M = 1- In turn, this alleviates the statistics require- 
ment (discussed later) by a factor of ~ 3. 

In Fig. 2, we display the survival probability (15) for 
the modern concordance cosmological parameters, i.e. 
fl M = 0.3, Q k = 0, fl A = 0.7, and h = 0.71 (top), 
and their allowed variations (bottom), respectively. It 
seems that variations of Om, Qa, and h, within their 
uncertainties [10, 11], amount to a ~ 5 % effect. 

As is apparent from Fig. 2 (top), using the narrow- 
width approximation (13) for the annihilation cross- 
section rather than taking into account the complete en- 
ergy dependence of the cross-section is justified within 
the overall 5 % error. 



B. Neutrino source emissivity distributions 

In the previous subsection, we have found that, for a 
given source emissivity distribution C V/3 of neutrinos of 
flavor P, the neutrino flux of flavor a to be observed at 
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FIG. 2: The survival probability P Vi [Ei,z) of a cosmic neu- 
trino i^i, injected at redshift z with energy as a function of 
the energy at Earth, E = Ei/(l + z), in units of the resonance 
energy E T J: B — m% /2 m Vi . 

Top: The narrow-width approximation (15), for z — 2 (dot- 
ted), z = 5 (short-dashed), z = 20 (long-dashed), and 
standard cosmological parameters (Qm = 0.3, Qa = 0.7, 
h = 0.71), compared with the complete energy dependence 
from the annihilation cross-section of Ref. [3] (solid). 
Bottom: The survival probability for z — 2 and standard 
cosmological parameters (solid) compared with the most ex- 
treme variations allowed by up-to-date global fits: Qm = 0.20, 
Q A = 0.78, h = 0.81 (dashed) and Q M = 0.40, Q A = 0.61, 
h = 0.62 (dashed-dotted) [11]. 



We will focus mainly on sources which produce pions 
("hadronic" sources), be they astrophysical Zevatrons 
(bottom-up) or non-accelerator (top-down) ones. From 



the decay sequence ^ — > /i = 



+ 2u u 



(nei- 



ther we, nor experiments, will distinguish neutrinos from 
anti- neutrinos), the flavor ratios of the source emissivities 
are predicted to be 



£„ : £„„ :£.. 1:2:0. 



(18) 



In this case, one finds - exploiting the fact that phe- 
nomenologically \U e 3\ 2 <§; 1 and {U^l ~ \U T z\ - that the 
fluxes at Earth are well approximated by (see Appendix 
and also Ref. [33]) 



1 

-in 



dz 1 



H{z) 3 



Ct ot (z,E(l + z)) (19) 



3 

£ 



\U aj \ 2 P Uj (E{l+z),z), 



where is the total, flavor-summed neutrino emissiv- 
ity at the source. In fact, as discussed in the Appendix, 
Eq. (19) holds whenever the source emissivities satisfy 
£„ fl +£„ T = 2£„ c . Thus, it holds also for a "democratic" 
flavor ratio of source emissivities, 



c ■ c ■ C — 1 • 1 ■ 1 

*^i/ e • *^v u • *-'l/ r X.X.I. 



(20) 



as might arise from the decay of topological defects not 
coupled directly to ordinary matter such as, for example, 
mirror-matter "necklaces" [50]. 

It is unlikely that neutrino observatories will be able to 
distinguish neutrino flavors at extreme high-energy [51]. 
Phenomenologically then, we are mainly interested in the 
sum over all neutrino flavors a. With such a sum, uni- 
tarity completely removes the dependence on the leptonic 
mixing matrix from Eq. (19), leaving just 



1 

4-7T 



dz 1 
H(z)3 



Ct ot (z,E(l + z)) 



Earth is predicted to be (cf. Eqs. (6) - (9)) 

oo 

5 /j£j x (17) 



]T \U aj \ 2 P Vj (E (1 + z), z) \U tij \ 2 (z, E(l + z)) , 

0,3 

where C Vg (z 1 E(l + z)) is the number of neutrinos of fla- 
vor [3 and energy Ei — E(l + z) emitted per co-moving 
volume, per unit time and unit energy, at a redshift "dis- 
tance" z (cf. Eq. (7)). In this subsection and the one 
that follows, we evaluate this expression further for some 
benchmark source emissivities. 



3 

x^2P Vi (E(l + z),z). (21) 
i=i 

Another simplification of Eq. (19) occurs if the neu- 
trino masses are quasi-degenerate, which in fact is re- 
alized in Nature if the lightest neutrino has a mass 
m Vl > Am| tm , say, m Ul > 0.1 eV (cf. Fig. 3): 

m Vl : m„ 2 : m„ 3 ~ 1 : 1 : 1 =^ P Ul ~ P V2 ~ P U3 . (22) 

With quasi mass-degeneracy, Eq. (19) simplifies to 

F„ a {E)~ (23) 

OO 
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FIG. 3: Allowed ranges for the neutrino masses as a func- 
tion of the lightest neutrino mass m Ul , in the normal (top) 
and inverted (bottom) three-neutrino scheme (adapted from 
Refs. [52, 53]). 



for each a — e, fj,, t. 

In the following, we assume that the EHECz/ sources, 
which build up the diffuse source emissivity, have iden- 
tical injection spectra J V0 ■ Consequently, the z and E$ 
dependences of the source emissivity distribution factor- 
ize, 

C V0 (z,E i )= V [z)J V0 {E i ). (24) 

The "activity" T](z) is the number of sources per co- 
moving volume and per unit of time, at the redshift 
"time" z, while the injection spectra J v „(Ei) are the 
number of neutrinos vp emitted by a single source per 
unit of energy. The z-dependence of the activity accounts 
for any evolution of the co-moving number density and/or 
of the common luminosity of the sources. 

We employ several paramctcrizations of rj(z), which 
allow us to study a variety of possible EHECV origins - 



ranging from astrophysical accelerator sources such as 
gamma ray bursts (GRB's) and active galactic nuclei 
(AGN's) which turned on at z ~ few, to non-accelerator 
sources such as topological defects which have been de- 
caying all the way back to z very large. 

We start with a parametrization of the activity moti- 
vated by astrophysics [54], 

7 ' SFR(Z) = 110 (l + z)-^a(l + z)»3 < ^ 

with r]o = rj{z — 0) being the activity in the today's 
epoch. With the values a = 0.005(0.0001), m = 3.3(4.0), 
and ri2 — 3.0(3.0), the parametrization fits the star for- 
mation rate (SFR) history derived from the blue and UV 
luminosity density, in line with the extreme ranges of 
optical/UV measurements without [55] (with [56]) dust 
extinction corrections. We will refer to these two cases 
as conservative and generous SFR, respectively. The 
star formation rate is believed to map out the earliest 
structures sufficiently bound by gravity to undergo fu- 
sion. As such, they may map out the history of AGN 
and GRB evolution, two potential sources of the EHEC 
rays (EHECR's) [54]. 

The parametrization (25) provides a peak activity at 



/ 72,1 \ ™l+™2 

1 + z poak = — , (26) 

\an 2 J 

which occurs at z pca k = 1.4 for the conservative SFR and 
at Zp 0a k = 2.9 for the generous one. For n 2 = 0, a = it 
reduces to a simple power-law, J](z) = tjq (1 + z) ni . Such 
a choice, with n\ ~ 4 to fit the low- 2 SFR, is often used 
in the literature. Without a cutoff in z, however, the 
power-law provides an extreme evolution history, as the 
activity increases indefinitely. 

This brings us to a further parametrization, namely a 
simple power-law ansatz, but with cutoffs z m j n and z max 
to exclude the existence of nearby and early-time sources, 
respectively: 

?7 P ow(z) = 770 (1 + z) n 9(z - z mhl ) 6(z max - z) . (27) 

This ansatz has the advantage of leading to easily trac- 
table analytic expressions, while still reproducing, for 
n = 4, the generous SFR case, as long as z max < z pca k. 
In addition to describing the evolution of GRB's and 
AGN's up to z ~ 2, the power-law with n ~ 1 2 and 
z mm ^> 1 also characterizes the activity expected from 
non-accelerator sources. For example, topological defect 
sources [57] are characterized by n = 3/2 [93] and z max 
arbitrarily large [94]. Throughout, we will take z min = 
default value [95]. 
Turning to the source injection spectra J V/3 (Ei), we 
only need to specify them in the energy decade around 
the respective resonance energies. For this energy decade, 
we again assume a power-law behavior [96] , 

Jvp{Ei) = j V0 E i " 6(Ei — E m i n ) 9(E max — E{) . (28) 
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As with the activity power-law, this spectral power-law 
parametrization facilitates the analytic study of different 
scenarios for the origin of EHECVs. 

In most of the source models, the neutrinos originate 
from pion decays, the latter either being produced in in- 
elastic pp or p7 interactions (astrophysical sources), or, 
alternatively, arising in the fragmentation of QCD jets in 
the decays of superheavy particles (top-down sources). 
Although the resulting neutrino spectra 3 Vfi can be fairly 
well calculated for given injection spectra of protons or 
pions, the details are tedious and not worth the effort 
for our present purposes, and so we retain the sim- 
ple parametrization (28) instead. The power-law should 
mimic the various neutrino injection spectra reasonably 
well, with the spectral index in the range a ~ 1 -j- 2. 
Throughout, we will take E m i n — as a default value. 

For astrophysical accelerators, one expects that -E max 
is a fraction (<~ 5 %) of the maximum proton energy, 
£ pmM . In the case of shock acceleration, E pmax is deter- 
mined by the requirement that the gyromagnetic radius 
of the protons in the ambient magnetic field B be less 
than the spatial extension R of the accelerating source. 
The result is 

Ef^l ~ 10 21 eV (R/kpc) (B/mG) . (29) 

Even higher energies are possible in proposed non- 
shock mechanisms, such as unipolar induction, acceler- 
ation in strong electromagnetic waves in plasmas (wake- 
fields) [58] , or by magnetic recombination in the vicinity 
of massive black holes [59]. 

For top-down scenarios, E maK can be much larger, ba- 
sically bounded only by the huge mass of the decaying 
particle or defect. These huge masses are thought to 
reflect the energy-scale of an underlying phase transi- 
tion. Popular examples include Mqut ~ 10 16 GeV from 
grand-unification, and M w i mpz in a ~ lO 11 ^" 13 GeV from 
the end of inflation. Neutrinos from the decay of these 
super-massive objects carry an energy about one order of 
magnitude less than the mass. 

In general, there may be several classes of sources - 
such as GRB's, AGN's, EHEC protons scattering in- 
elastically off the CMB ( "cosmogenic" neutrinos), and 
topological defects - which build up the total emis- 
sivity distribution. Predictions from a variety of such 
source classes are shown in Fig. 1 (bottom): from jets 
of AGN's [60] , from ordinary topological defects (Mx = 
10 14 GeV) [61], from hidden-sector topological defects 
(M x = 4 x 10 14 GeV) [62], and from the Z-burst sce- 
nario normalized to fit the highest energy cosmic ray 
anomaly [24] with different assumptions about the uni- 
versal radio background. Also shown are the upper and 
lower bounds [46] (short-dashed-shaded) and one exam- 
ple [61] (short-dashed) of the cosmogenic neutrino flux. 
Finally, the cascade limit [63] from Ref. [64] (dotted- 
shaded) on transparent neutrino sources (discussed later) 
is shown. It applies to all scenarios where neutrinos origi- 
nate from pion decays or even from electroweak jets [65]. 
These neutrinos are accompanied by photons and elec- 



trons which cascade down in energy during their prop- 
agation through the universe. The cascade limit arises 
from the requirement that the associated diffuse gamma- 
ray fluxes should not exceed measurements [97]. 

It is straightforward to generalize Eq. (24) to a sum of 
source classes. In reality, it is more prudent to wish for 
even one class of sources at the energies of interest here, 
i()22h-23 eV, and so we continue to work within the single 
source-class hypothesis. 

Another possibility, in principle, is to collect events 
from a small number of point sources, possibly just one. 
This presents the advantage that the spectra are not 
smeared by the integration over the redshift distribu- 
tion. Since neutrinos are not deflected in flight, point- 
source selection is possible. However, we will find that 
the event numbers required for statistical significance is 
on the order of 100 or more. This large number presents 
a luminosity challenge for a small number of very bright 
sources [2]. In the end, it does not matter much for our 
purposes whether the EEC^ flux is diffuse of granular, 
since the isotropy of the relic background ensures nearly 
the same absorption shapes in the spectrum of either 
source- type. 



C. Case studies 

Let us start our analyses with the quasi-degenerate 
neutrino mass-spectrum case, whose implications are 
summarized in Eq. (22). In view of the expectation that 
EHEC^ neutrino fluxes are rapidly falling with energy 
(cf. Fig. 1), this case has the best prospects of observ- 
ability, since it corresponds to the largest neutrino masses 
(cf. Fig. 3), and it brings all three resonance energies 
to a common lowest value (1). As discussed above (cf. 
Eq. (23)), we expect in this case an identical absorption 
dip in the neutrino flux of each flavor to be observed at 
Earth, if the flavor ratios at the sources were hadron- 
like (18) or democratic (20). 

In Fig. 4, we show the predicted flux F„ a for a hadronic 
or flavor democratic source emissivity characterized by 
a power-law activity (27) with redshift evolution index 
n, and power-law injection spectrum (28) with index 
a. This flux is normalized to the predicted flux for no 
absorption, F"° abs (obtained by replacing the survival 
probability in Eq. (23) by unity). A nice feature follow- 
ing from the two power-laws, source activity and source 
emissivity, is that the normalized spectrum depends on 
n and a only through the combination n — a: the source 
evolution and the energy fall-off compensate each other 
in a simple way. Three particular n — a combinations are 
shown in Fig. 4 for each fixed z max : n — a = (upper 
curves), n — a = 2 (middle curves), and n — a = 4 (lower 
curves). 

Figure 4 illustrates some general features. Viewed as 
a function of decreasing energy, the absorption dip starts 
abruptly at the resonance energy E™ s . This initial depth 
is determined by today's annihilation probability, «3% 
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FIG. 4: Predicted flux (23) of neutrinos u a of flavor a — e,fj,,r 
at Earth, for a source emissivity characterized by a power-law 
activity (27) and a power-law injection spectrum (28), for the 
case of quasi-degenerate neutrino masses (22), normalized to 
the predicted flux for no absorption. E / Ef B scales as the 
degenerate mass ra„. Three values of z m ax are shown: 2 
(dotted), 5 (short-dashed), and 20 (long-dashed). For each 
choice of 2 max , three choices of n — a are shown: (upper), 2 
(middle), and 4 (lower). The corresponding solid lines show 
the same quantity evaluated with the complete energy depen- 
dence of the annihilation cross-section from Ref. [3] arising 
from the finite Z-width, instead of exploiting the zero-width 
approximation (13). For all curves, £ max > E r J^ (1 + z mal ) is 
assumed. 



according to (16), and reflects the absorptions occur- 
ring in nearby (z ~ 0) space. The dip extends in en- 
ergy down to E™ s / (\ + z max ), with this minimum en- 
ergy being the redshifted value of the resonant energy 
for annihilations occurring at cosmic distance [z ~ z max ). 
The overall depth of the dip increases dramatically with 
n — a, thereby showing a strong preference for source 
evolution and/or flat energy-spectra. With increasing 
z ~ z max = 2, 5, 20, the absorption depths are roughly 5, 
8, and 15 % for n — a = 0; 7, 18, and 55 % for n — a = 2; 
and 10, 27, and 77 % for n — a = 4, respectively. We 
observe in Fig. 4, as we did in Fig. 2 (top), that the 
replacement of the finite Z-width with the 5-function ap- 
proximation (13) is well justified for our purposes. From 
now on, we will exploit this simplification. 

Next, we consider some non-degenerate neutrino mass 
scenarios, and source activities other than power-law. In 
the four panels of Figs. 5-8, we show the prediction (21) 
f° r Fv a /J2 a abs for a neutrino spectrum which is 
quasi-degenerate, i.e. lightest mass m Vx = 0.4 eV (upper 
panel), normal hierarchical, with m Vl = 0.01 eV (2nd 
panel) and m Vx — 0.002 eV (3rd panel), and inverted 
hierarchical with m Vl = 0.002 eV (bottom panel). Fig- 
ures 5 and 6 display results for z max = 10, a — 2 (solid) 
and (dashed), and the conservative and generous star- 
formation activities, respectively. Figure 7 shows results 



for a power-law activity, n — a = 0, and z max = 10 (short- 
dashed) and 20 (solid), mimicking a topological defect 
source scenario. Fig. 8 displays results for a power-law 
activity, n — a = 2, and z max = 2, 5, 10 from upper to 
lower curves, corresponding to bottom-up acceleration. 
As expected, the depths of the dips increase markedly 
with z max , lower spectral index a (28), and increased 
source evolution (27), the latter represented either by 
the shift from "conservative" to "generous" SFR as in 
Figs. 5 and 6, or by the increase in n — a as in Figs. 7 
and 8. 

The 2nd panels in the figures illustrate that in the case 
of a normal hierarchy and an intermediate lightest neu- 
trino mass of m Vl = 0.01 eV, corresponding to central 
values to„ 2 = 0.013 eV and m„ 3 = 0.053 eV in Fig. 3 
(top), the absorption dips from v\ and v-i cannot be re- 
solved. The two dips appear as one dip, which is un- 
derstandably about twice as deep as the single dip from 

The two lower panels present results for a very small, 
lightest neutrino mass of m„, = 0.002 eV. In the nor- 
mal hierarchical case (3rd panel), three absorption dips 

- associated with the masses m Vl = 0.002 eV, and, from 
the central values of Fig. 3 (top), m„ 2 = 0.0085 eV and 
m„ 3 = 0.052 eV - are clearly visible. These 3rd pan- 
els present relic neutrino mass absorption spectroscopy 
at its best! For an inverted hierarchy, on the other 
hand, one finds, as illustrated in the bottom panels, 
that the absorption dips from U2 and 1/3 cannot be re- 
solved individually. With m Vx = 0.002 eV, the cen- 
tral values for the nearly-degenerate heavy masses are 

= 0.052 eV, from Fig. 3 (bottom). 
The 2 meV mass is a critical one in the sense that 
the relic neutrino v\ is indeed non-relativistic for z up 
to z max = 2, in which case Eq. (12) is applicable. For 
even smaller masses (or larger z max ), the relic v\ density 
becomes relativistic and the absorption dip is known to 
wash out [1]. Thus, the two lower panels in Figs. 5-8 
remain applicable for any m Vl < 0.002 eV, if one simply 
removes the contribution of the highest-energy dip, at 

- 10 24 eV. 



III. EXPERIMENTAL PROSPECTS 

In the last paragraph we have seen that, depending 
on the source activity, injection spectrum, and neutrino 
masses, we may expect absorption dips of (10 -7- 20) % 
depth, located at (0.1 -j- 0.5) £v s , with a width of about 
an order of magnitude in energy. Is there any hope of 
discovery of these absorption dips in the next decade or 
beyond? The answer to this question depends critically 
on the magnitude of the EHECV flux at the resonance en- 
ergies realized in Nature. In the following we will study 
a few benchmark flux scenarios (§ III A) and discuss fur- 
ther experimental (§ III B) issues concerning the clean 
reconstruction of absorption dips. 
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FIG. 5: Predicted flux of neutrinos summed over flavors at FIG. 6: Same as Fig. 5, but with the generous SFR activ- 
Earth (21), normalized to the predicted flux for no absorption, ity (25). 
for a conservative SFR activity (25) to z ma x = 10, injection- 
spectrum indices a — 2 (solid) and a — (dashed) (28), and 
neutrino spectra which are quasi-degenerate (top), normal- 
hierarchical (2nd and 3rd panels), and inverted-hierarchical 
(bottom panel). For all curves it is assumed that .E max > 
Elf{l + z mBix ). 
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FIG. 7: As in Fig. 5, but with a power-law activity (27) and 
n — a = 0, with 2 m i n = and z max = 10 (short-dashed), 
20 (solid), mimicking a topological defect source scenario. For 
all curves it is assumed that £ maj( > E r °* (1 + z max ). 



FIG. 8: Same is Fig. 7, but with n — a = 2 fixed, and 
^max = 2,5,10 (from upper to lower curves), corresponding 
to a bottom-up acceleration source scenario. 
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TABLE I: Expected number of neutrinos (plus anti-neutrinos) 
to be detected in upcoming EHECV observatories with ener- 
gies in the indicated intervals until the indicated year, for two 
different EHECV flux scenarios - one saturating the current 
observational upper bound and one saturating the cascade 
limit (cf. Fig. 1 (bottom)). 







energy decade 


10 21 - 22 eV 


1Q 22-23 e y 


10 23 - 24 eV 


year 


2008 


2013 


2008 


2013 


2008 


2013 


observ. limit 


240 


700 


30 


90 


2 


5 


cascade limit 


13 


40 


3 


10 


1 


2 



A. Benchmark flux scenarios 

1. Most optimistic scenario: hidden sources 

The most favorable case for relic neutrino absorption 
spectroscopy is the one in which the neutrino fluxes satu- 
rate the current observational upper bounds. It is known 
quantitatively (cf. the points marked "Z-burst" from 
Ref. [24] in Fig. 1 (bottom)) that in this most favorable 
flux scenario, the secondary fluxes of protons (and pho- 
tons) from hadronic Z-decay are remarkably of just the 
right order of magnitude to explain the highest energy 
cosmic rays above the GZK energy by Z-bursts [20, 22]. 

As can be seen from Fig. 1 (bottom) and as summa- 
rized in Table I, the upcoming EHECV observatories ex- 
pect to see in this case in 2013 about 230 neutrinos (plus 
anti-neutrinos) per flavor a — e, [i, r in the energy inter- 
val from 10 21 eV to 10 22 eV. The total number of neutrino 
events, then, is N ~ 700 neutrinos of all flavors, which 
implies a 1-sigma fluctuation of y/N ~ 26. For a 3-sigma 
evidence for an absorption dip in this energy interval, an 
absorption depth of 3 vN /N ~ 11 % is required; for a 5- 
sigma discovery, a depth of 5 \J~N /N ~ 19 % is required. 

As a comparison with our case studies in § II C re- 
veals, these depth requirements are achievable as long 
as the dip is maximized via a quasi-degenerate neutrino 
mass. The quasi-degenerate condition is m Ul > 0.1 eV, 
which in turn implies that E™ s <4 x 10 22 eV. For exam- 
ple, power-law source emissivities produce such depths 
for n — a > and/or z max > 10 (cf. Figs. 4, 7 (top), and 8 
(top)). Therefore, the realization of the Z-burst mecha- 
nism for the already observed highest energy cosmic ray 
events implies not only a discovery of the corresponding 
EHEO flux by 2008 (cf. Table I), but also the discov- 
ery of solid evidence for the associated absorption dip as 
soon as 2013. 

What kind of source could produce such an opti- 
mistic EHECz^ flux, of the order of the observational 
limit at 10 21 eV to 10 22 eV? Bottom-up astrophysical 
Zevatrons may in principle produce such a flux, with 
the problem that they must accelerate protons to ener- 




10 18 10 19 10 ao 10 ai 10 aa 10 23 



Ea E [eV] 

FIG. 9: Predicted neutrino flux at Earth, summed over all 
flavors, from a power-like source emissivity, with n = 1.5, 
•Zmax = oo, a = 1.5, and £ max = 4 x 10 22 eV. This 
flux mimics one from hidden-sector topological defects with 
M x = 4 x 10 14 GeV (cf. Fig. 1 (bottom)) and is also suf- 
ficient to explain the EHECR's above -Egzk via the Z-burst 
mechanism. Curves are without (dotted) and with relic neu- 
trino absorption. Assumed neutrino masses are degenerate at 
m Vl = 0.2 eV (dashed) and m Vl = 0.4 eV (solid). The error 
bars indicate the statistical accuracy achievable per energy 
decade by the year 2013, for a flux which saturates today's 
observational bound from Fig. 1 (bottom). 



gies E p max>10 23 eV [24, 25]. Much more favorable in 
this respect are top-down sources, involving physics be- 
yond the Standard Model (SM) which naturally produce 
extremely-high energetic particles. With either bottom- 
up or top-down, the sources must be "hidden" [98], in 
order to avoid the cascade limit (cf. Fig. 1 (bottom)). 
Hidden sources are, by definition, sources from which 
neither nucleons nor photons (with energies > 100 MeV) 
escape. 

One popular hidden top-down example is a topological 
defect, here taken with mass Mx ~4x 10 14 GeV, which 
couples to SM particles only indirectly through a non-SM 
sector [50]. In this case of topological-defect origin, we 
have n — a ~ and z max ^> 1, and so Fig. 7 applies; 
the dip presented there is ~ 15 % for the most inter- 
esting case of quasi-degenerate neutrinos, m Vx > 0.1 eV. 
In Fig. 9, we show the significant wiggle from quasi- 
degenerate neutrinos in the otherwise power-law spec- 
trum expected for the hidden topological-defect neutrino 
flux as in Fig. 1 (bottom). By design, this flux scratches 
the observational upper bound. Such a flux is also suffi- 
cient to explain the EHECR's via the Z-burst mechanism. 
The indicated error bars show the statistical significance 
that is expected with planned and proposed experiments 
by the year 2013 (cf. Fig. 1 (bottom)). This example 
illustrates concretely our previous claim that substantial 
evidence for a relic neutrino absorption dip and the ex- 
istence of the C^B is achievable within the foreseeable 
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FIG. 10: As in Fig. 9, but with _E max = 10 24 eV, to mimic 
topological defects with Mx = 10 16 GeV (cf. Fig. 1 (bottom)). 
The assumed neutrino spectra are: (i) quasi-degenerate, 
m vi = 0.4 eV (solid), (ii) normal hierarchical, m vi = 0.01 eV 
(long-dashed) and m vi = 0.002 eV (short-dashed), and (hi) 
inverted hierarchical, m vi = 0.002 eV (dashed-dotted). Here 
the error bars indicate the statistical accuracy achievable per 
energy decade by the year 2013, for a flux which saturates 
today's cascade limit from Fig. 1 (bottom). 



future, if the Z-burst mechanism for the EHECR's is re- 
alized in Nature. 

For a hierarchical neutrino spectrum, on the other 
hand, with m Vl < 0.04 eV, the lowest resonant-energy is 
~ 10 23 eV. Even if the EHEC^ flux saturates the cur- 
rent observational limit, the expected event numbers in 
the 10 22 + 23 eV energy interval are apparently too small 
(cf. Table I) to allow a discovery of absorption dips in 
the next decade: with 90 expected events in this energy 
range by 2013, one needs absorption depths of ~ 53 % 
(~ 32 %) for a 5-sigma discovery (3-sigma evidence). As 
can be gleaned from Fig. 7, such large dips are not ex- 
pected. 



2. Less optimistic scenario: transparent sources 

Let us turn our attention now to the non-optimized 
EHECz; fluxes. From a neutrino flux that saturates 
the cascade limit (cf. Fig. 1 (bottom)), we may ex- 
pect, in the 10 2lH ~ 22 eV energy bin, just 40 events by 
the year 2013 (cf. Table I). This yields a 3-sigma ev- 
idence for an absorption dip only if it has a depth of 
~ 48 % [99]. Dips of this depth we have seen in § II C 
only for quasi-degenerate neutrinos and extreme parame- 
ter choices for the source activity. For example, for power 
indices n — a > 2, as might happen with highly- evolving 
bottom-up sources, one needs z max > 20 in addition (cf. 
Figs. 4 and 8 (top)). The latter is easily achieved in top- 
down sources (z max = oo), however their restriction to 
n — a ~ tends to decrease the depth of the dip (cf. 



FIG. 11: As in Fig. 10, but with a generous SFR activity (25), 
Zmax = 20, injection spectrum index a — 2, and _E max = 
10 24 eV, to mimic astrophysical sources with -E pma x = 2 x 
10 16 GeV. 



Figs. 4 and 7 (top)). This pessimistic outlook is how- 
ever ameliorated when one realizes that an increase in 
statistics by a factor of 10 reduces the required absorp- 
tion depth by a factor of ~ 3, to ~ 15 % for a 3-sigma 
evidence, and to ~ 25 % for a 5-sigma discovery. Such 
an increase in statistics could be achieved, for example, 
by undertaking more ANITA flights, by extending the 
EUSO flight time, or by developing the OWL or SalSA 
experiments. If the ~ 13 neutrinos in the lO 21 "^ 22 eV in- 
terval (cf. Table I) , expected from a flux which saturates 
the cascade limit, are measured by the pre-2008 exper- 
iments, then such extensions of the upcoming EHECV 
observatories would be warranted. 

What classes of sources could deliver EHECV fluxes 
which saturate the cascade limit? It has recently been 
pointed out that cosmogenic neutrino fluxes from photo- 
pion production by cosmic protons on the CMB can reach 
this limit [46, 61]. However, the limit is scratched typi- 
cally just above -Eqzk, at around 10 20 eV and not higher 
(cf. Fig. 1 (bottom)). Other possible sources include 
bottom-up Zevatrons with a large E pmax > 10 23 eV, and 
topological defects with M x > 10 14 GeV. 

As examples, we show in Figs. 10 and 11, for various 
neutrino mass spectra, the predicted modulations in the 
otherwise power-law spectra expected for a topological 
defect and an astrophysical neutrino flux, respectively, 
both saturating the cascade limit. As expected, we see 
that a further increase in statistics by a factor of ~ 10 -j- 
100, to reduce the error bars by a factor of ~ 3 ~ 10, will 
be mandatory to establish absorption dips. Moreover, a 
quasi-degenerate neutrino mass spectrum is required. 
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B. Further experimental issues 

To temper optimism, it is fair to mention three further 
issues that mitigate an experiment's ability to cleanly 
observe an absorption dip. These are (i) the resolution 
with which an experiment can reconstruct the initial neu- 
trino energy from the visible event energy, (ii) the flavor- 
dependent nature of energy reconstruction, and (iii) the 
possible confusion of a deviation of the (assumed) power- 
law spectrum due to absorption with other possible ori- 
gins of deviation. We discuss each of these in turn. 

Proposed EHECV detectors will measure shower en- 
ergy from ground-based scintillator or water, or atmo- 
spheric nitrogen fluorescence, or radio signals in ice. 
From EHECR physics, energy reconstruction for hadron- 
initiated showers is known to be about 25 %. Neutrino 
observatories expect to reconstruct shower energies with 
about the same error. This will smear the dip somewhat, 
but our analysis is based upon events per decade in en- 
ergy, so the 25 % smearing should be tolerable. 

However, there is a subtlety connected with the various 
showers that result from different neutrino flavors [100]. 
Consider first neutral current (NC) events, identical for 
all neutrino flavors. On average, the final state neu- 
trino carries away 80 % of the incident energy, leav- 
ing 20 % in the detected jet. The NC cross-section is 
<~ 45% that of the charged-current (CC) cross-section. 
Next consider CC interactions. The and v T CC events 
produce charged leptons, generally unobservable, which 
carry away 80 % of the energy. Thus, the and v T 
CC events leave the same 20 % energy deposition in the 
shower as do the NC interactions. For all these events, 
the absorption dip in the shower events will appear lower 
in energy than the true dip by a factor ~ 0.2. This is 
easily corrected. However, there remains some smearing 
due to the event-by-event variance of the energy trans- 
fer about the mean. We do not include the effect of this 
variance in the present work. 

In contrast to the above event classes, a v e CC event 
will produce a hadronic jet plus an electron which creates 
an electromagnetic jet, and so the interaction deposits 
100 % of the incident energy into the combined shower. 

In summary of this flavor discussion, NC, and and 
y T CC scattering constitute about 77 % of the events. 
The showers from these events contain about 20 % of the 
incident neutrino energy The other 23 % of events are 
due to v e CC, with all of the incident energy observed. 
Thus, the observed events will look like a superposition 
of two fluxes with relative weights 23 % and 77 %, the 
latter displaced downward in energy by a factor of five in 
the mean, but with fluctuations. 

Finally we come to the confidence issue, whether an 
observed feature near the end of the EHEC^ spectrum 
can be claimed to be an absorption dip. For example, the 
spectral end-point may have structure simply due to dif- 
ferences among the individual contributing sources. To 
better ensure that a dip feature is observed, the contin- 
uation of the spectrum above the dip region should be 



observed. This requires more events, at higher energies 
still. Looking again at Table I, one sees even more reason 
to doubt the observation of an absorption feature, except 
for relatively low resonance-energies. In turn, this means 
that quasi-degenerate neutrinos, with mass not much be- 
low present cosmological bounds, are required if the relic 
neutrinos are to be detected via an absorption dip. 



IV. SUMMARY AND CONCLUSIONS 

Relic neutrino absorption spectroscopy via the obser- 
vation of absorption dips in the EHEC^ flux may be fea- 
sible. The Z-dips, if observed, are rich in particle and as- 
trophysical information. The position of the high-energy 
edge of each dip is fixed by the neutrino mass. The depth 
and shape are determined by the density of the CjaB, by 
the usual cosmological parameters, by the activity and 
injection spectrum of the EHEC^ sources, and by the 
neutrino mass pattern. Moreover, as we have seen, the 
width of the dip region reflects the spectrum and evolu- 
tionary history of the neutrino source(s). Thus, were we 
so fortunate as to resolve a dip in some detail, the in- 
formation to be mined truly belongs to interdisciplinary 
particle- astrophysics. 

However, the statistics at the dip is entirely determined 
by the magnitude of the EHECf flux. Moreover, the en- 
ergy positions of the dips - and, to some extent, their 
depths - critically depend on the neutrino mass spec- 
trum. 

Large event samples, N » 100, beyond 10 21 eV are 
needed to reveal Z-dips with statistical significance. To 
get these event numbers within the next decade or so, an 
EHECi^ flux at least as large as the present cascade limit 
is required. Almost certainly, even higher fluxes must 
be invoked. They could be generated by hidden sources, 
which are opaque to nucleons and high-energy photons, 
thereby evading the cascade limit. Such sources are likely 
indicative of new physics. Moreover, the high flux re- 
quirement implies a Z-burst contribution to the EHECR 
events beyond Eqzk within an order of magnitude of 
the present AGASA rate. The Auger project should de- 
tect these large EHEC^ and EHECR fluxes within a few 
years, and EUSO should easily measure them within the 
decade. 

As inferred from our numerous figures, a quasi- 
degenerate neutrino mass spectrum with m Vl > 0.1 eV 
seems to be required to produce a detectably-deep ab- 
sorption dip at a possibly-accessible resonant energy. 
Such a spectrum is testable in several ways [9] . Neutrino- 
less double beta decay experiments (assuming neutrinos 
are Majorana) such as CUORE [74] and NEMO-3 [75], 
and the KATRIN tritium decay experiment [76] can be 
expected to show positive results already in the upcom- 
ing decade. The transfer functions relating CMB fluc- 
tuations to today's large-scale matter distributions can 
be expected to show positive contributions from massive 
neutrinos [77]. It seems sensible to say that, if pre-2008 
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experiments do not see any EHECV flux in the 10 21 ^" 23 eV 
region, then, in the context of the concordance cosmolog- 
ical model, absorption dips won't be observed within the 
next decade or two. If such is the case, then Nature will 
have overlooked a wonderful opportunity to produce di- 
rect evidence for the C^B [101]. 

In summary, the presently planned neutrino detectors 
open up a window of opportunity for relic neutrino ab- 
sorption spectroscopy. The next decade will be really 
exciting and decisive in this respect. 



APPENDIX: IMPLEMENTING NEUTRINO 
FLAVOR PHYSICS 

At production in the cosmic sources, let the ratios of 
the neutrino flavors be written as v e : : v T = w e : 
: w T , with Y^, w = 1- A convenient description of the 



flavor mixture is given by the density matrix 

P{t = 0) = ^wp\ Vp )( V p\. (A.l) 



The density matrix is properly normalized to Tr(p(0) = 
1, and so describes the ensemble-averaged, single neu- 
trino. 

The relic neutrinos, with thermal energies ~ 3kT u ~ 
0.5 meV in the today's epoch, long ago decohered into 
their mass-eigenstates [21]. Therefore, rewriting the den- 
sity matrix in the mass-basis allows us to simply include 
the losses due to resonant absorption. Making use of the 
mixing matrix notation \vp) = Upj \ or equivalcntly, 
Ufjj = {vj\vfj) — , the forward-propagated den- 

sity matrix in the mass basis, with resonant absorption, 
is: 



p(t) = X 



E 



]T wp ]T Ufa U* 0k e+ iA < 4 e~ r -^ 4 \vi) {v k \ . 

j,k 



(A.2) 



Here, Tj = enj Oj is the annihilation rate, and the factor ergy losses due to redshifting; however, it is straightfor- 
e -5 1 corr ectly accounts for absorption at the amplitude ward to incorporate redshifting in the final expressions, 
level [102]. As written, this formula does not include en- At Earth, the probability to detect flavor a is then, 



P Va detected = (v a | p(t) K> = w S UfijUtkUakKj e~ iA "^ * * . (A.3) 

I 



This probability expression naturally divides into a sum from which oscillations arise. After a bit of algebra, the 
diagonal in the mass, and a sum of interfering mass terms result is: 

I 



^detected = Wf, ^ \U aj \ 2 \Upjf e" 1 * * (A.4) 

3 

x - x - r j+ r , 

+ 2 2^ W 0Z^ e 2 

j<k 

r 



^{u N u; k u ak u* ak ) cos 



2E 



( U 03 U 0k U <*k U aj) sin 



2E 



Setting the annihilation rates to zero returns the usual formula for neutrino oscillations (cf., e.g., Ref. [53]). 
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Although the neutrino wave-function at extreme high- 
energies may remain coherent over a cosmic distance, the 
oscillating terms are effectively averaged away in any re- 
alistic observation. Let us look first at the coherence of 
the neutrino wave-function. The spread in the neutrino 
mass-eigenstates after travel through a distance D results 
from the difference in the group velocities (3 = SE/Sp: 

D Am 2 

Spread = Ct 5 f3 = ^ E 2 ( A ' 5 ) 

= 6 x 10- 20 ^ (JL^j Am 2 _ 3 E^ cm, 

where Am 2 _ 3 = Am 2 /T(r 3 eV 2 , E 22 = E/10 22 eV, and 
the Hubble distance D H = cH^ 1 = 4.2 (0.7 /h) Gpc char- 
acterizes the typical cosmic distance. Decoherence of the 
wave packet results when L> sp read exceeds the natural 
length of the wavepacket, call it ct$, the latter being 
determined by conditions at the source [103]. The natu- 
ral length may be c times the production time [80] (e.g. 
CTj± ~ 3 m), the mfp between interactions in a dense 
source [80], or the spatial uncertainty in the location of 
the production point in the source [81]. The decoherence 
length is obtained by setting Asprcad in Eq. (A. 5) equal 
to ct$, and solving for D. The result, as a fraction of 
the Hubble size, is 

gdecohere = (h \ / T^X 1Q 20 (A fi) 

D H \0.7 J V3m/ Am 2 3 v ' 

Clearly, decoherence does not occur in our Universe for 
a 10 22 eV neutrino. 

Now we turn to the averaging effects of measurement. 
To observe the oscillating terms requires measurement of 
the phase <f> = Am 2 t/2E to better than 2ir, or equiv- 
alcntly, knowledge of S(D/E) to better than Air/Am 2 . 
This in turn requires knowledge of 5D to better than 
Aosc = At: E/ Am 2 , and knowledge of SE to better than 
E (X osc /D). With the oscillation length being so short 
<~ E22/ ' Am 2 _ 3 kpc compared to cosmic scales, there is 
no hope to observe the oscillating term. 

From here on, we may use for the detection probability 
just the averaged value of Eq. (A. 4) , i.e. just the first 
term 

P Va detected - l^l 2 e ~ T] ' Yj W P l^"! 2 ' ( AJ ) 

3 

Finally, taking into account energy losses due to redshift, 
Eq. (A.7) is easily generalized to Eq. (9) in the main text 
for the differential propagation function of neutrinos, and 
to Eq. (17) for the neutrino flux to be observed at Earth. 

As simple as Eq. (A.7) is, further simplifications are 
possible. For example, if w e : : w T = 1 : 1 : 1, as 
might arise from "democratic" neutrino emission from a 
topological defect, then by unitarity of the mixing ma- 
trix, wp \ U@j | 2 — 1/3 independently of j. It is also 



known [33] that the same result obtains for the case 
w e : Wfj_ : w T = 1 : 2 : as results from charged pion 
decay to neutrinos, in the limit of <-> v T interchange- 
symmetry. The latter is exact if the atmospheric mixing 
angle #32 is maximal (45°) and ^St(U e s) = 0. Oscillation 
data [5] show that these two conditions are satisfied (or 
nearly so). Since any linear combination of the two sets 
{wp} just discussed must also give the same conclusion, 
we arrive at a mini-Theorem: 

whenever w e = 1/3 and <-» v T interchange-symmetry 
is (nearly) valid, 

then ^2gW/3 \U{3j\ 2 — 1/3 independently of j. 
Because w e — 1/3 includes two popular neutrino- 
production cases, charged pion decay and democratic 
emission, we will assume that w e = 1/3 holds and work 
with the the very compact result 

Pv a detected = 3 X/ \Uaj\ 2 j * . (A.8) 
3 

This expression leads immediately to the main text's 
Eq. (19) for the neutrino flux at Earth, after again taking 
into account the energy loss due to redshifting. 

Cosmic-neutrino detectors are unlikely to be capable 
of neutrino flavor identification (sec Ref. [51] for a dis- 
cussion). Accordingly, it is tempting to also sum over a, 
the flavor arriving at Earth, to arrive at the especially 
simple result 

P ^ y = E P ^ detected = \ ]T e~ T > * , (A.9) 
a j 

leading to Eq. (21). Caution is warranted, however, as 
the interactions of the different flavors deposit different 
mean energies into the detector. For example, the short 
mean-free-path (mfp) of a high-energy electron, and the 
long mfp's of high-energy muons and taus, imply that 
in charged-current interactions, j^'s and i^ T 's will leave 
much less visible energy in the detector than will v e 's. 
This issue is discussed further in the main text in § III. 
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spin-states per flavor). 
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[99] If the number of events are small, one should apply Pois- 
son statistics, as given in Ref. [72]. 

[100] We take our numbers from the extrapolated cross- 
sections in [73]. 

[101] Even if Z-dips cannot be measured, it may still be pos- 
sible to infer the C^B from Z-burst data. The statis- 
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as predicted by the concordance model. The parameters 
of the concordance cosmology may give way to some- 
thing new. All of this encourages an open mind. If large 
structures did form early, then more neutrino sources 
and strong source evolution may be the reality. With 
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